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There is increased interest from the research commu- 
nity and chnicians to implement closed-loop stimulation 
strategies in neurobionic devices. That is, to adjust sti- 
mulation levels dynamically based on the responses of 
neural tissue in real time. To adjust electrical stimula- 
tion in a closed-loop bionic device, a model-based con- 
troller design can be implemented. Here, we collect 
experimental data from retina slices and use data-driven 
technique to model neural dynamics. Our motivation 
comes from visual prostheses. 

In vitro experiments were conducted on NZ white rab- 
bit retina tissue. Electrical stimulation of the retinal gang- 
lion cells consisted of a train of pulses whose amplitudes 
had a white Gaussian distribution with zero mean and 
standard deviation of luA. In this pulse train, there were 
a total of 5000 biphasic pulses with 100|is/phase dura- 
tion, equal phase amplitude, and no interphase delay. 
Frequency pulse trains of 25-lOOOHz were used. 

Experimental data were used to characterize the 
response properties of neurons to the electrical stimula- 
tion. The reverse-correlation approach, widely used to 
predict neural response to light stimulation, was adapted 
to predict neural response of ganglion cells to electrical 
stimulation. In contrast to the traditional reverse-correla- 
tion schemes, the model proposed here incorporated the 
history of the response of a neuron. To emphasize this, 
we call it the "spike history" model. To estimate the 
responses of cells to electrical stimulation, we used a 
novel pseudo-random stimulation. To validate the fitness 
of the model, we performed statistical analysis of the 
simulated spike trains. In particular, we compared the 
values of the coefficient of variation of the interspike 
interval for the experimentally recorded spike train, for 
the simulated spike train with the history model, and for 
the simulated spike train using the model without the 
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response kernel. To compare how well the simulated 
spike train approximated experimentally recorded spikes, 
we compared the penalty terms for the spike history 
model and for the model without the response kernel. 
The penalty term was calculated based on the time differ- 
ence between each simulated spike and the closest spike 
in time in the experimentally recorded train. 

Cells responded to all stimulation frequencies. A phe- 
nomena of clusters of spikes followed by periods of sup- 
pression was observed in raw spike trains above 200 Hz 
stimulation. The phenomenon of cluster-suppression 
with white-noise stimulation may be indicative of mem- 
ory in the system; i.e., the response depends not only on 
the current stimulus but also on the responses within a 
time window preceding the current time. To confirm 
this, we calculated the auto-correlation function of the 
recorded spike trains for different frequencies of stimula- 
tion. This showed that there was memory in the system 
for stimulation at frequencies higher than lOOHz. 
Robustness of the model parameters was confirmed by 
repeated stimulation of the same cell at the same fre- 
quency (up to 20 repetitions). Statistical comparison of 
the model with history and the model without history 
confirmed that the model with history approximates 
experimental data better. 

The proposed model can be used to design a stimula- 
tion strategy to control neural excitation. 
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